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Abstract 

We revisit the question of quantum dynamics of electrons on the off-diagonal Fibonacci 
tight-binding model. We find that typical dynamical quantities, such as the probability 
of an electron to remain in its original position as a function of time, display log-periodic 
oscillations on top of the leading-order power-law decay. These periodic oscillations with the 
logarithm of time are similar to the oscillations that are known to exist with the logarithm 
of temperature in the specific heat of Fibonacci electrons, yet they offer new possibilities for 
the experimental observation of this unique phenomenon. 

1 Fibonacci Electrons 

The Fibonacci sequence of Long (L) and Short (S) intervals on the 1-dimensional line — generated 
by the simple substitution rules L — )• LS and S — )• L — is a favorite textbook model for 
demonstrating the peculiar nature of electrons in quasicrystals [IHS]. The wave functions of 
Fibonacci electrons are neither extended nor exponentially-localized, but rather decay alge- 
braically; the spectrum of energies is neither absolutely-continuous nor discrete, but rather 
singular-continuous, like a Cantor set; and the quantum dynamics is anomalous. In recent years 
we have studied how these three electronic properties change as the dimension of the Fibonacci 
quasicrystal increases to two and three [IHH] , by constructing square and cubic versions of the 
Fibonacci quasicrystal [9j. 

The 1-dimensional off-diagonal Fibonacci tight-binding model is constructed by associating 
a unit hopping amplitude between sites connected by a Long interval, and a hopping amplitude 
T > 1 between sites connected by a Short interval, while assuming equal on-site energies that 
are taken to be zero. The resulting tight-binding Schrodinger equation, on an Fjv-site model, is 
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given by 

T,-+i^(i + l) + T,V'(j-l) = i?V(j), j = l...FN, (1) 

where ip{j) is the value of an electronic eigenfunction on site j, E is the corresponding eigenvalue, 
-Fat is the N*-^ Fibonacci number, and the hopping amplitudes Tj are equal to 1 or T according to 
the Fibonacci sequence {Tj} = {1, T, 1, 1, T, 1, T, 1, 1, T, 1, 1, T, 1, T, 1, 1, T, 1, T, . . .}, as described 
above. The diagonal Fibonacci tight-binding model is constructed by taking all the hopping 
amplitudes to be equal to 1, and taking the on-site energies to have two values, Vl and V5, 
arranged according to the Fibonacci sequence. These models have been studied extensively ever 
since the initial interest in the behavior of electrons in quasiperiodic potentials [lOHlSj . and 
continue to offer mathematical challenges to this day [14]. 

Various 2-dimensional extensions of the Fibonacci model were introduced soon thereafter [151 - 
I19j . and strongly promoted recently [9] as models for quasicrystals without 'forbidden' symme- 
tries [20p21j. In our studies, we have shown that whereas Fibonacci electrons in one dimension 
always behave as described above for any T > 1, in two dimensions, and even more so in three, 
there is crossover — as the strength T of the quasiperiodicity approaches 1 — to a regime in which 
Fibonacci electrons behave more and more like electrons do in periodic crystals, particularly 
in the sense that their energy spectra develop continuous intervals [IHS]. These results were 
recently explained in a rigorous manner by Damanik and Gorodetzki [22]. More surprisingly, 
our studies of Fibonacci electrons have led us to new results also in the simple 1-dimensional 
case. We have examined dynamical properties, such as the probability of an electron to remain 
in its original position as a function of time. The power-law decay of this quantity is commonly 
used for analyzing the dynamics. Surprisingly, we have observed log-periodic oscillations on top 
of the power-law decay, implying the existence of an imaginary correction to the exponent. We 
wish to describe this new observation here. 



2 Quantum Dynamics of Electronic Wave-Packets 

We consider the dynamics of electronic wave-packets, or states |n), that are initially localized at a 
single lattice site n of the 1-dimensional Fibonacci quasicrystal, denoting their amplitude on site 
m at time t by (pniiT^, t)- Thus, at time t = the wave-packet is given by (pni^, 0) = (m|n) = 6mn, 
where 6mn is the Kronecker delta. At any later time t > the wave-packet is given by e*^*|n), 
which is simply the n*'^ column of the matrix representation of the time evolution operator e*''^*, 
where H is the off-diagonal matrix representation of the Hamiltonian in Eq. ([1]), and we take 
n = l. Thus, (l)n{m,t) = {m\e'^^n) = {e'^*)mn- 

We characterize the dynamics of wave packets by monitoring two typical quantities: (a) The 
Survival Probability of the n^^ wave-packet, defined as the probability of finding the electron at 
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its initial position at time t, is given by 
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Sn{t) = \Mn,t)\' = (e^^Onn ; (2) 

and (b) The Inverse Participation Ratio of the n*^ wave-packet, which measures the spatial 
extent of the wave-packet, is given by 

m = = 5: |(e-*)m. ' , (3) 



(Em \M'm,t)\ 



where the last equality holds because the wave-packets are normalized. Both of these quantities 
are often used to examine the dynamics of wave-packets. In particular, the manner in which 
they decay as the wave-packets spread with time is associated with different regimes of the 
quantum dynamics. The actual calculation is performed numerically by applying the discrete 
time-evolution operator e^^^^ successively to the initial conditions. The time-step At is kept 
sufficiently small to satisfy the Nyquist criterion by noting that the largest eigenvalue of H is 
bounded by 1 -|- T, thus requiring At to be smaller than 1/(2 + 2T). 

If a function F{t) decays asymptotically with some power law F{t) ~ t~^ , then the exponent 
(3 can be found by 

t^co In t 

but in general the exponent j3 is not guaranteed to exist. The bounds I3± on /3, which always 
exist, are given by 

(5+ = - limsup — /3_ = - liminf — 5 

and if /3+ = /3_ then /3 exists. Furthermore, if an exponent < /3 < 1 characterizes the 
power-law decay of a function F{t)^ then the same exponent also characterizes the decay of the 
time-averaged function 

{F), = - f F{t')dt'. (6) 
t Jo 

We use the definition of Eq. ([5]) to study the long-time asymptotic values of the exponents /3s 
related to the survival probability of a wave-packet and /3/ related to the inverse participation 
ratio of the wave-packet. In what follows we are mostly concerned with the early-time behavior 
of these exponents. 

Damanik et al. |14^I23) recently used the second moment of the position operator for the 
diagonal tight-binding Hamiltonian, to show that far from the periodic limit, or very close to it, 
the dynamics of wave-packets is independent of the initial site. However, since we are interested 
in intermediate values of T as well — equivalent in the diagonal model to intermediate values of the 
difference |Vl — between the different on-site energies — we expect to find different dynamical 
behavior depending on the choice of the initial site. Thus, we typically examine the maximal, 
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Figure 1: The time-averaged maximal survival probability for a 233-site Id approximant, with 
periodic boundary conditions, calculated for different values of T. The asymptotic behavior of 
the slopes allows to extract the diffusion exponent f3g^°-^ as a function of T, which is displayed 
in Fig. [21 

the minimal, and the site-averaged survival probabilities and inverse participation ratios, all 
of which display similar qualitative behavior. In Ref. [8] we study the different exponents in 
detail, in one, two, and three dimensions. Here we concentrate on the Id results for the maximal 
survival probability exponent 13^°'^ . 

3 Asymptotic Behavior of the Maximal Survival Probability 

We look at the wave-packet whose time-averaged survival probability is maximal. Fig. [T] shows 
the time-averaged maximal survival probability {Smax)t on a log-log scale for a 233-site Id ap- 
proximant, with periodic boundary conditions, for different values of T. The exponents extracted 
from the slopes of the curves in Fig. [1] for four orders of approximants are shown in Fig. [21 as 
functions of T. 

Convergence of the exponent /3^°^ is evident for values of T > 2 even for relatively small 
approximants. Convergence is not obtained for smaller values of T that approach the periodic 
limit T ^ 1, where the dynamics is expected to become ballistic. In this limit the wave- 
packets quickly spread out and the finite size of the approximant has a stronger influence on 
the dynamics. Despite this limitation, detailed studies of the slopes, for increasing orders of 
approximants and for varying time-ranges clearly show that the extracted exponents converge, 
as the order of the approximant increases, and longer averaging times are possible. We note 
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Figure 2: The exponent /3g^"^ as a function of T for increasing orders of approximants. Conver- 
gence is improved for decreasing values of T as the order of the approximant is increased. 

that Similar curves are obtained for the minimal and the site- averaged survival probabilities, as 
well as for the minimal, maximal, and site-averaged inverse participation ratios [8j. 

4 Log-Periodic Oscillations and the Fourier Transform of the 
Density of States 

A closer inspection of the temporal decay of the survival probabilities in Fig. [1] reveals a small 
oscillating behavior on top of the overall leading power-law. In order to better expose this be- 
havior, we divide out the leading asymptotic power-law that was found earlier, and obtain 
the curves shown in Fig. [3l Similar results are observed in the study of the inverse participa- 
tion ratio. The curves clearly exhibit log-periodic oscillations — oscillation that are periodic in 
log(t) — around the mean value of 1, especially for the larger values of T. These oscillations 
contain a basic frequency u that seems to decrease with increasing T, as well as a sequence 
of higher-frequency oscillations with decreasing amplitudes that seem to develop with time. 
The fundamental oscillations can be described empirically by a temporal decay with a complex 
exponent of the form 




(7) 



which after dividing out the leading power-law term become 



fit) 



cx 1 + at" 



1 + ae 



,±iu] log(t) 



(8) 
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Figure 3: Log-periodic oscillations in the maximal survival probability, after dividing out the 
leading power-law decay, for an 89-site approximant with 5 different values of T. In addition 
to the basic periodic oscillation one observes the emergence of higher-frequency oscillations as 
time advances. 
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yielding oscillations of amplitude a around 1 that are periodic in log(t) with frequency uj. 

Log-periodic oscillations are not as uncommon in physics as one may think. They appear in 
critical phenomena |24fl26] . where renormalization-group calculations yield power-law behavior 
with complex exponents near phase transitions; they also appear quite generally in problems 
related to random-walks or classical diffusion on self-similar systems or fractals, where the fractal 
dimension turns out to be complex [27ti3T] . As such, they naturally occur in quasiperiodic 
systems like the Fibonacci quasicrystal, where the spectrum of energies is multifractal. This has 
been observed in quasiperiodic Ising models [32], and in specific heat studies of the tight-binding 
Fibonacci model [331437] . 

In specific heat studies, where log-periodic oscillations appear in temperature, the object of 
calculation is the partition function 



Z(/3t) 



k 



(9) 



where /3t is the inverse temperature; whereas in problems of diffusion the object of calcu- 
lation is the so-called heat kernel of the diffusion equation 



Z{t) = ^ e 



(10) 



In both cases, when expressed in integral form, these are Laplace transforms of the density 
of states. To see the relation to the present problem of quantum dynamics, we expand the 
expression for the survival probability ([2]) in eigenstates of the Hamiltonian, 



Sn{t) 



n e 



\-^k){'^k\n) 



(11) 



Thus, what we are calculating is the magnitude-squared of the Fourier transform of the density 
of states, weighted by the overlap of each eigenstate with the initial state |n). 

Equation (|lip is easily calculated in the periodic limit T — ?> 1 of our model, where the 
eigenstates are Bloch function, with \ipk{'n)f' = 1/N independent of n, and the eigenvalues form 
a single band with dispersion = cos(/c). This yields 



ST=l{t) oc 



dke 



it cos k 



OC J2(t), 



(12) 



where Jo{t) is the zeroth order Bessel function of the first kind, which is known to decay asymp- 
totically as This is consistent with the expected ballistic dynamics of Bloch electrons in a 
periodic crystal. The time average of Eq. (|12p is approximately the behavior calculated numer- 
ically for r = 1 on a 233-site approximant, and shown in Fig. [TJ Note that the fine wiggles that 
are barely seen at very early times are associated with the first few zeros of the Bessel function. 
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and are not log-periodic oscillations. Log-periodic oscillations appear for T > 1 as a result of 
the multifractal nature of the Fibonacci spectrum, for similar reasons as in the calculation of 
the partition function ([9]) |361l37j. A more detailed analysis of Eq. (jlip for T > 1 is required 
to quantitatively characterize the spectral properties of the log-periodic oscillations, which we 
have discovered here by numerical means. 

Log-periodic oscillations in the specific heat of Fibonacci quasicrystals near zero temperature 
are very difficult to observe experimentally. The discovery of log-periodic oscillations in the 
quantum dynamics of wave-packets in Fibonacci quasicrystals should open new possibilities for 
the actual experimental observation of this unique phenomenon. A particular realization could 
be in optical experiments that allow one to observe the dynamics of wave-packets within 1- 
dimensional [38], as well as 2-dimensional [39], photonic quasicrystals. Thus, we hope that our 
numerical observations here will stimulate further studies, analytical and experimental alike. 
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